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U.  ABSTRACT 

»*  v 

'A  computer  program  has  been  written  to  calculate  the  electric  field  re¬ 
quired  to  initiate  microwave  breakdown  in  a  partially  ionized  gas.  The  break¬ 
down  field  is  found  by  calculating  iouization  rates  and  diffusion  loss  rates  and 
then  solving  for  the  value  of  the  electric  field  at  which  these  rates  balance. 

the  preionization  case,  a  modified  diffusion  coefficient  is  required  to  ac¬ 
count  for  charge-separation  effects.  The  program  reported  here  uses  the 
method  of  MacDonald  for  calculating  ionization  and  diffusion  loss  rates  and 
extends  that  method  to  the  preionization  case  by  using  a  multiplicative  factor 
to  model  the  transition  to  ambipolar  diffusion.  Results  are  graphed  for  neon 
and  neon-argon  mixtures  using  various  values  of 'the  characteristic  diffusion 
length  and  signal  frequency  over  a  range  of  gas  pressures. 
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Microwave  Breakdown  Calculations  That  Include 
the  Effects  of  Preionization  in  Neon  an  i 

Neon -Argon  Mixtures 


L  INTRODUCTION 

To  calculate  the  electric  field  requir-d  to  initiate  microwave  breakdown  in  a 
partially  ionized  gas,  a  computer  program  was  written  and  an  expression  derived 
for  the  rate  of  electron  loss  due  to  diffusion  in  the  free  and  ambipolar  modes.  The 
breakdown  calculations  reported  here  are  based  on  a  simplified  form  of  this  diffu¬ 
sion  expression  and  therefore  include  the  effects  of  charge  separation  resulting 
from  preionization.  Ionization  by  electron-metastable  collisions  and  other  mech¬ 
anisms  that  may  become  important  in  dealing  with  preionization  are  not  treated 
here. 

1.1  The  Breakdown  Criterion 

Microwave  breakdown  occurs  in  a  gas  when  the  rate  of  electron  production 
equals  or  infinitesimally  exceeds  the  rate  of  electron  loss.  The  electric  field 
accelerates  the  electrons  which  subsequently,  through  collision  with  neutral  gas 
atoms,  randomize  this  ordered  oscillatory  energy.  Since  energy  is  coupled  into 
the  electron  constituent,  an  electron  avalanche  will  occur  if  a  sufficient  number 
of  electrons  accumulate  a  random  kinetic  energy  greater  than  the  ionization  poten¬ 
tial  of  the  neutra.  gas  atoms.  In  this  case,  ionizing  collisions  of  electrons  with 


(Received  for  publication  5  September  1973) 
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neutrals  will  occur,  electron  multiplication  will  result  when  the  ionization  rate 
exceeds  the  loss  rate,  and  progressively  more  electrons  will  be  present  to  be  ac¬ 
celerated  by  the  electric  field.  In  microwave  breakdown,  which  is  the  sudden  on¬ 
set  of  a  high-frequency  discharge,  the  ambient  electron  density  is  generally  low 
enough  for  the  electron  loss  to  be  by  diffusion.  Where  production  equals  loss, 

.  t 

the  breakdown  condition  can  be  stated  as  an  electron  continuity  equation: 

=  ni/.  -  V  •  T  =  0  ,  (1) 

where  u  is  the  ionization  frequency  and  T  is  the  electron  flux  to  the  wa^ls.  If 
the  diffusion  coefficient  D  is  defined  by  setting  the  flux  equal  to  the  gradient  of  a 
diffusion  potential,  as 

?  =  -  V  (nD),  (2) 


then  the  continuity  equation  can  be  rewritten  as 

|f  =  nv.  +  V1 2 3  (nD)  =  0.  (3) 

This  equation  is  separable,  and  eigenvalues  can  be  f'.und  for  the  separation  con¬ 
stant  by  considering  boundary  conditions.  It  is  convenient  to  denote  these  eigen¬ 
value  solutions  by  A,  a  characteristic  diffusion  lengt  -  whose  value  for  various 

2 

geometries  is  determined  by  the  dimensions  of  the  dis  charge  walls. 

3 

Using  A  yields  the  breakdown  condition  (analogous:  to  the  Townsend  criterion  ): 


(4) 


Next,  the  breakdown  field  is  solved  for  by  calculating  V.  and  D  for  various  elec¬ 
tric  fields  and  iterating  onto  the  electric  field  that  satisfies  (4). 

1.2  The  Ste.dy-Statc  Cue 

In  a  low-density,  quiescent  microwave  discharge,  electron  diffusion  losses 
are  reduced  by  the  space-charge  fields  set  up  by  an  excess  positive  charge  of 


1.  Brown,  S.C.  (1966)  Basic  Data  of  Plasma  Physics,  M.l.T.  Press,  pp.  251,  252. 

2.  MacDonald,  A.  D.  (1966)  Microwave  Breakdown  in  Gases,  Wiley,  pp.  27-29. 

3.  MacDonald,  A.  D.  (1968)  Microwave  Breakdown  in  Gases,  Wiley,  pp.  69-72. 
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ions  'left  behind'  by  their  more  rapidly  diffusing  electrons.  Electron  losses 
occurring  in  the  transition  from  free  electron  diffusion  to  restricted  ambipolar 

diffusion  are  described  by  a  space-charge-controlled  diffusion  coefficient  D  , 

4  8 

defined  in  the  steady  state  by 


Ds  = 


(5) 


The  approximation  to  D  given  by  Allis  and  Rose  is  the  constant- ratio 

s 

approximation 


D  *  D 
s  r 


(6) 


where  r  is  the  ratio  of  positive-ion  concentration  to  electron  concentration,  and 
p.  is  the  mobility  of  species  i . 

In  the  space-charge-controlled  case,  where  r  =  1 ,  the  usual  arrLipolar- 
diffusion  coefficient  D  is  obtained  by  reducing  D  ,  as  follows . 

3  *  I* 


/u  D_  +  M.dA 

°r  =  \  **+  +  <»-  j5  °a‘ 


In  the  low-density  free-diffusion  case,  where  r  =  D_/D+  ,  the  exp-ession  for  Dr 
reduces  to  just  D_ ,  In  the  intermediate  range,  a  machine  solution*  has  shown 
that  the  approximation  underestimates  Dr  <  Dg  by  a  factor  of  no  more  than  3 . 
Manipulating  the  expression  for  Dr  yields  the  more  convenient  form: 


D_  =  D 


(7) 


To  program  this  requires  a  suitable  function  for  r ;  it  must  approach  D  /D  for 

2  2  2  2  -  + 

Xp  »  A  and  unity  for  XD  «  A  ,  where  XD  is  the  Debye  length.  This  is  satis¬ 
fied  by  the  following,  which  is  chosen  for  convenience. 


4.  Allis,  W.  P. ,  and  Rose,  D,  J.  (1954)  The  transition  from  free  to  ambipolar 
diffusion,  Phys,  Rev.  93(No,  l):84-'"3. 
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In  Figure  1,  where  r  is  plot*“d  as  a  function  of  ( A/X^j  it  is  observed  to  be  pro¬ 
portional  to  (A/Xj^  in  the  transition,  which  is  proportional  to  the  electron  den¬ 
sity.4  Substituting  Eq.  (8)  into  (7)  for  Dr  yields 


(9) 


Figure  1.  An  Analytical  Form  for  r  is  Necessary  for  Programing. 
Equation  (8)  has  been  chosen  to  model  r  and,  as  shown  here,  has  the 
correct  limits,  as  well  as  a  linear  region  proportional  to  electron  density. 
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Thic  expression  seems  less  awkward  when  reduced  to  a  familiar  form  for  a 
specific  case.  Consider  the  isothermal  case,  where  D _/|i_  =  D+/|i+.  Now 


(10) 


■* 

If  D+/D_  in  the  denominator  is  neglected  as  small  with  respect  to  unity,  then  for 
the  high-density  case  (where  «  A2)  the  term  Dr  reduces  to  the  ambipolar  dif¬ 
fusion  coefficient: 


D  =  2D  s  D.  ..  ,  . 

r  +  isothermal 


In  this  case  of  microwave  breakdown,  where  the  electron  temperature  far  ex¬ 
ceeds  the  ion  temperature,  such  a  restrictive  isothermal  approximation  cannot  be 
used.  The  ambipolar  coefficient  is  then  obtained  by  considering  the  high-density 
case  in  Eq.  (9),  from  which 


Kp.  +  M-P A 

\  +  r 


(id’ 


This  can  be  reduced  by  neglecting  a-  sm.iil  with  respect  to  p_  .  Using  the 
Einstein  relation  then  yields 


(12) 


again,  a  familiar  form. 
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1.3  The  Preionization  Caae 

As  used  here,  the  term  preionization  denotes  the  state  of  a  gas  that  has  been 
partially  ionized  by  an  energy  source  independent  of,  and  consequently  acting 
prior  to  application  of,  the  microwave  breakdown  signal.  It  is  assumed  that  the 
density  of  charged  particles  in  the  initial  ionization  that  exists  as  an  afterglow, 
for  example,  or  as  a  weakly  ionized  discharge  is  sufficient  to  require  considera¬ 
tion  of  charge-separation  effects  and  subsequently  a  diffusion  coefficient  that  is  a 
function  of  electron  density. 

5.  Oskam,  H.  J.  (1958)  Microwave  investigation  of  disintegrating  gaseous  dis¬ 
charge  plasmas.  Philips  Res,  Rpts  13:367. 
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This  breakdown  crite'.*ion  (Eq.  (4)  with  a  modified  D]  predicts  that  breakdown 
will  occur  for  the  preionization  case,  in  the  presence  of  an  electric  field  lower 
than  the  field  in  the  ambient  case.  The  expression  used  for  modeling  the  preioni¬ 
zation  case  represents  a  constant-ratio  approximation  diffusion  coefficient. 


2.  PROGRAM  FUNCTIONS  AND  FUNCTION-TESTING 

g 

Bandel  and  MacDonald  have  shown  experimentally  that  the  field  required  for 
microwave  breakdown  in  neon  is  reduced  by  preionization.  In  my  calculations  of 
the  breakdown  fields  in  neon  and  neon-argon  mixtures,  I  use  a  multiplicative 
factor  to  allow  for  the  change  in  diffusion  coefficient  caused  by  preionization. 

7 

The  breakdown  field  for  neon  has  been  calculated  by  MacDonald  and  Betts, 
who  used  a  Boltzmann  equation  solution  for  the  electron  velocity  distribution  at 
both  high  and  low  pressures.  The  breakdown  field  for  neon-argon  mixtures  has 

O 

been  similarly  calculated  by  MacDonald  and  Matthews. 

The  computer  program  numerically  calculates  values  for  the  distribution 
functions  in  the  analytic  expressions  as  they  appear  after  integration  over  electron 
velocities,  and  yields  the  ionization  and  diffusion  loss  rates.  For  neon  at  low 
pressures,  modified  Bessel  functions  of  imaginary  arguments  of  the  first  and 
third  kinds  are  required.  Foi  the  other  cases,  the  confluent  hypergeometric 
functions  are  required.  [The  second  solution  given  by  MacDonald6 7 8 9  also  requires 
the  r  and  t|>  functions  for  its  evaluation.  ] 

The  modified  Bessel  functions  were  initially  calculated  from  ordinary  Bessel 
functions  of  the  first  and  second  kind  by  using  imaginary  arguments.  It  was  later 
found  that  direct  calculations  using  series  and  asymptotic  representations  are 
more  efficient. 

A  simple  program  was  written  Co  test  the  functions.  Each  function  subroutine 

_4  v 

was  called  for  various  values  of  its  argument  (ranging  from  10  to  10"),  and  the 
value  of  the  function  wa3  assigned  to  an  array.  After  all  the  values  of  the  argu¬ 
ment  had  been  computed,  the  arrays  were  printed  out  and  checked  against  tables 
for  the  various  functions.  The  exponential  integral  function  and  the  modified  Bes¬ 
sel  function  of  the  first  kind  were  accurate  to  at  least  the  eighth  place.  The 


6.  Bandel,  H.  W. ,  and  MacDonald,  A,  D.  (1969)  Effect  of  preionization  on  micro- 

wave  breakdown  in  neon,  J.  Appl,  Phys.  40(No.  ll):4390-4394. 

7.  MacDonald,  A.  D. ,  and  Betts,  D.  D.  (1952)  High-frequency  gas  discharge 

breakdown  in  neon,  Canad,  J.  Phys.  30:565-576. 

8.  MacDonald,  A.  D. ,  and  Matthews,  J.  H.  (1955)  High-frequency  ionization  co¬ 

efficients  in  neon-argon  mixtures,  Phys,  Rev.  98(No.  4):  1070- 1073. 
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modified  Bessel  function  of  the  third  kind  was  also  accurate  to 'eight  places,  ex¬ 
cept  in  the  range  of  arguments  around  15 .  where  the  algorithm  changed  from  the 
series  representation  to  the  asymptotic  form.  Here  the  values  were  accurate  to 
only  four  places.  The  confluent  hypergeometric  functions  were  checked  against 

g 

the  limited  range  of  values  given  by  MacDonald  to  the  fourth  place  and  were 
found  to  agree  exactly.  For  values  outside  this  range,  the  asymptotic  represen¬ 
tations  were  used.  These  agreed  with  NBS  tables  to  the  third  significant  figure, 
but  it  was  found  that  the  arguments  of  the  functions  used  in  the  breakdown  field 
programs  were  never  in  the  asymptotic  range  except  possibly  in  the  first  few 
iterations  for  a  particular  case. 


3.  BREAKDOWN  CALCULATIONS 
3.1  ‘Fake  Position’  Iteration 

The  field  required  to  initiate  breakdown  was  calculated  by  using  a  version  of 
the  'false  position'  iteration  technique.10  In  this  technique,  two  initial  values  of 
the  argument-in  this  case,  the  electric  field-are  assumed.  They  are  then  tested 
to  find  out  whether  the  solution  lies  between  them.  If  it  does,  a  difference  equa¬ 
tion  calculates  a  'closer'guess,  which  is  in  turn  tested,  and  so  on.  The  method 
hinges  on  writing  a  function  that  has  a  solution  when  the  function  is  equal  to  zero. 
The  product  of  the  function  evaluated  at  the  two  trial  (false  position)  values  is 
negative  if  the  solution  lies  between  these  two  values,  positive  if  the  solution  lies 
outside  the  bounds  of  these  trial  values. 

Recall  our  breakdown  criterion: 


This  equation  can  be  rearranged  into  such  a  function,  as  follows: 


((4)] 


(13) 


A  specific  f(x)  was  written  for  each  of  the  four  cases:  for  neon,  at  high  and 
low  pressures;  and  for  neon-argon,  at  high  and  low  pressures.  The  iteration 
program  was  the  same  for  all  four  cases. 


9.  MacDonald,  A.  D.  (1949)  Properties  of  the  confluent  hypergeometric  function, 
J.  Math.  Phys.  28:183-191. 

10.  Conte,  S.D.  (1965)  Elementary  Numerical  Analysis,  McGraw-Hill,  pp.40-41. 


11 


Th ;  value  of  the  gas  pressure  was  incremented  after  each  value  of  the  elec¬ 
tric  field  had  been  iterated  onto  for  the  various  valucc  of  electron  density  at  that 
gas  pressure.  After  the  electric  field  had  been  calculated  for  breakdown  without 

O 

preionizaiion  (in  the  program,  n  =  1. 0cm*'5),  suitable  values  of  electron  density 
were  used  to  span  the  region  of  transition  from  free  diffusion  to  ambipolar 
diffusion. 


Returning  to  Eq.  (91  and  setting 

A  ‘  T\  (,2'2'!c.kT-) 

H  D.  +  M.D  /  22l  \ 

B  =  T^rnrw+ 


we  can  rewrite  (9)  as: 

Dr  *  D-(&Sl)-  <16> 

The  magnitude  of  factors  A  and  B  determines  the-,  electron  density  at  which 
transition  occurs.  The  ratio  A/B  represents  the  reduction  in  diffusion  by  space 
charge  as  expression  (16)  approaches  A/B  for  large  n.  The  ratio  of  A  to  B  cSn 
be  estimated  from  Eqs.  (14)  and  (15),  which  yield: 


5,  a  a*  i  a-  iv 

D_  "  B  "  M+  +  I*. 

Neglecting  p+  as  small  with  respect  to  leads  to 

4-5.5:  „ 

B  fl_  +  D_  ' 

Since  the  ratio  of  the  diffusion  coefficient  to  the  mobility  is  equal  to  the  average 
energy,  then 

(*+/)*_  D_/M_  <  u_> 

dtd;=  <u+>  •  “ 
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In  the  laboratory  case  of  microwave  breakdown,  the  electron  average  energy 
is  much  greater  than  the  ion  average  energy.  Thus,  neglecting  D+/D_  as  small 
compared  with  M+^M_  lea.1s  to: 


In  the  Maxwellian  equilibrium  case,  using  the  Einstein  relation  gives  the 

ratio  of  electron-to-ion  temperatures  rather  than  the  average  energies  of  (19). 

For  purposes  of  the  program  described  here,  the  electron  mobility  was  esti- 

1 1 

mated  by  the  following  linear  approximation  to  the  data  of  Anderson 
v_  =  106  E/p, 
which  yielded 


-  -E 


p  V-sec  * 


For  the  ion  mobility,  the  value  4  crn  /V-sec  was  selected  since  experimental 
results12  give  values  for  small  E/p  slightly  larger  than  4  ,  decreasing  with  E/p. 
The  term  A  was  estimated  by  assuming  the  following: 

D_  "  ft_T_  ' 


.  T  =  300°K ; 

V-sec  + 

6  2 

-r£— —  ,  T  =3000°K, 
p  V-sec 


Including  the  constants  for  neon  yielded 


A  =  pA2(2. 8  X  10" 12). 


11.  Anderson,  J.  M.  (1964)  Hall  effect  and  electron  drift,  Phys.  Fluids  7:1517. 

12.  Chanin,  L.  M. ,  and  Biondi,  M.A.  (1957)  Temperature- dependence  of  ion 

mobilities  in  helium,  neon,  and  argon,  Phys,  Rev.  106:473-479. 
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Based  on  (20),  the  term  B  was  estimated  as: 

A(2.5X105) 

B  =  -  . 

The  ion  mobility  for  the  neon-argon  mixture  was  assigned  a  value  of 
2 

1.  S  cm  /  V- sec,  on  the  assumption  that  practically  all  the  ions  were  argon.  Since 
electrons  colliding  with  neutrals  would  still  encounter  neon  neutrals  to  the  extent 
of  the  concentration  of  neon  (99  percent),  the  electron  mobility  of  neon  was  used. 
Thus, 

A  =  pA2(l.  12  X  10"12), 

B  _  A(6. 25  X  105) 


3.2  The  Electron  Energy  Distribution 

The  Boltzmann  equation  solutions  for  the  distribution  function  given  by 

7  8 

MacDonald  and  Betts  for  neon,  and  by  MacDonald  and  Matthews  for  a  neon- 

argon  mixture,  were  used  to  calculate  the  ionization  rate  and  diffusion  coefficient 

for  particular  values  of  an  applied  electric  field.  These  calculations  were  then 

used  in  the  breakdown  criterion  [Eq.(13)  with  (16)  substituted]  to  iterate  onto  a 

breakdown  field. 

The  Boltzmann  equation  is  solved  by  assuming  that  the  mean  free  path  is  con¬ 
stant  for  electrons  of  all  energies  and  that  all  electron  losses  are  by  diffusion  to 
the  walls  of  the  container.  The  characteristic  diffusion  length  A  for  a  cylinder  of 
radius  R  and  length  L  is  found  by  using: 


(22) 


A  spherical  harmonic  expansion  of  the  distribution  in  velocity  and  a  Fourier 
series  expansion  in  time  leads  to  the  following  energy  balance  equation: 


l2  f  _  E2l2  _d_  uu2 
3A2  3u  du  v2+u? 


df_ 
du  ' 


(23) 


where 


h*  is  the  efficiency  of  excitation, 
h.  is  the  efficiency  of  ionization, 


i4 


f  is  the  electron  velocity  distribution  function, 
m  is  the  electron  mass. 

M  is  the  ion  mass, 
u  is  mv^/2e .  in  units  of  voltage, 
t  is  the  electron  mean  free  path, 

E  is  the  rms  amplitude  of  the  microwave  field, 

0)  is  the  angular  frequency  of  the  incident  microwave  field, 
V  is  the  momentum  transfer  collision  frequency. 


2  2 

This  equation  is  then  solved  by  considering  two  cases-low  pressure:  v  «  w  ; 
2  2 

and  high  pressure:  u)  «  V  .  For  each  case,  an  expression  is  derived  for  f : 

(1)  for  energies  below  the  excitation  energy  u^ ,  ignoring  inelastic  collisions,  and 
12)  for  energies  above  ux ,  considering  inelastic  collisions.  These  expressions 
and  their  derivatives  are  then  matched  at  u  =  ux  to  give  the  constants  in  the 
solutions. 

Next,  the  distribution  function  expressions  are  used  to  calculate  the  ionization 
rate  for  neon  by 


w 


i 


z  Vdu 


(24) 


and  for  the  Penning  mixture  of  neon-argon  by 


(25) 


For  both  neon  and  neon-argon,  the  diffusion  coefficient  is  calculated  by 


nD 


8rr-te2 

3m2 


ufdu , 


(26) 


The  ratio  of  nv.  and  nD  is  then  used  in  Eq,  (13)  with  Eq.  (16)  p'*oviding  a  multi¬ 
plicative  factor  to  represent  the  charge-separation  effect.  This  gives: 


FCN  = 


nvi 

nD_ 


nv 


.2  =  nD  (TFA) 


1 

7 


(27) 


where  FCN  is  the  function  f(x)  of  Eq.  (13),  and  TFA  is  the  multiplicative  factor 
representing  the  transition  from  _free  to  ambipolar  diffusion. 

These  calculations  were  separated  into  four  cases  and,  as  previously  noted 
(Sec.  3.1),  a  program  written  for  each  (neon,  low  pressure:  neon,  high  pressure; 
neon-argon,  low  pressure;  neon-argon,  high  pressure). 
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3.3  COMMOtl 

The  iteration  technique  used  in  this  work  tests  for  convergence  by  comparing 
each  new  value  of  the  electric  field  with  the  value  from  the  preceding  iteration. 
The  value  of  the  electric  field  is  assumed  to  have  converged  to  at  least  one  sig¬ 
nificant  figure  less  than  the  difference  between  two  successive  iterations.  For 
example,  if  two  successive  iterations  show  a  difference  of  no  more  than  one- 
tenth  in  the  value  of  the  electric  field,  the  result  is  assumed  to  have  converged 
in  the  units  place. 

A  false  indication  of  convergence  may  result  if  the  initial  false  positions,  or 
bounds,  for  the  value  of  the  electric  field  are  too  far  apart.  The  function  fix)  is 
therefore  printed  out.  If  the  printed-out  value,  which  should  be  zero  at  the  solu¬ 
tion,  is  not  small  enough  (approaching  zero),  the  initial  guesses  (false  positions) 
are  edited. 

Values  of  the  arguments  for  the  confluent  hypergeometric  functions  are  listed 
to  see  the  range  of  arguments  used  and  determine  whether  or  not  the  asymptotic 
series  is  being  used.  These  are  listed  separately  from  the  breakdown  results  to 
facilitate  checking  values  of  fix)  for  convergence  and  analyzing  the  breakdown 
results. 

The  values  of  the  characteristic  diffusion  length  and  microwave  frequency  that 
were  chosen  were  those  of  my  experimental  setup  and  the  published  experimental 
results  of  MacDonald.  The  results  for  the  case  without  preionization  are  pre¬ 
sented  in  Figure  2,  together  with  the  experimental  results  cited.  Figures  3  to  5 
present  the  results  for  the  case  with  preionization,  where  the  diffusion  loss  rate 
is  modified  by  the  function  described  in  Sec.  1.  2,  and  the  reduction  in  breakdown 
field  is  seen  as  a  function  of  electron  density. 

The  program  was  written  with  an  interactive  feature,  as  follows.  At  the  be¬ 
ginning  of  the  program  a  line  was  printed  at  the  teletype  terminal  requesting  a 
value  of  characteristic  diffusion  length,  followed  by  a  line  requesting  a  value  of 
frequency.  Consequently,  any  other  value  of  these  parameters  could  readily  be 
entered  through  this  procedure.  For  convenience,  the  ratio  between  the  break¬ 
down  value  of  the  elect  c  field  with  a  given  level  of  preionization  and  that  without 
preionization  was  also  calculated  and  printed  out. 

4  DISCUSSION 
13 

Brown  has  briefly  summarized  the  subject  of  high-frequency  breakdown 
(Chap.  10).  In  describing  the  steady-state  microwave  discharge  (Chap.  14)  he 

13.  Brown,  S.  C.  (1966)  Introduction  t?  Electrical  Discharges  in  Gases,  Wiley, 
Chap.  10  and  p.  221^ 
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Figure  2.  Experimental  and  Calculated  Breakdown  Field 
as  a  Function  of  Gas  Pressure 


presents  experimental  curves  for  hydrogen  (p.  228)  that  show  the  balance  of  elec¬ 
tron  losses  and  gains  in  the  discharge  is  maintained  by  a  lower  electric  field  than 
that  required  to  initiate  the  discharge  (see  Figure  6).  The  difference  is  due  to 
ambipolar  diffusion  in  the  sustaining  case  as  contrasted  with  free  diffusion  for 
breakdown  without  preionization.  When  a  microwave  signal  is  propagating  in  the 
partially  ionized  gas  of  a  reentry  plasma  sheath,  therefore,  breakdown  can  be 
expected  at  a  lower  field  strength  than  in  the  case  of  an  un-ionized  signal  path. 
'Breakdown'  is  understood  to  mean  that  the  electron  density  s  increased  to  the 
critical  density. 

Among  laboratory  experiments  showing  that  breakdown  fields  are  reduced 
with  preionization  are  those  reported  by  Bethke  and  Reuss.  These  workers 
have  observed  that  a  transient  microwave  breakdown  phenomenon  initiated  by  a 
glow  discharge  occurs  at  field  strengths  much  lower  than  the  field  strength 

14.  Bethke,  G.  W. ,  and  Ruess,  A.D.  (1969)  Microwave-induced  plasma  shield 
propagation  in  rare  gases,  Phys.  Fluids  12:822. 


17 


BREAKDOWN  VOLTAGE  (V/cm)  BREAKDOWN  VOLTAGE  (V/cm)  BREAKDOWN  VOLTAGE  (V/cm) 


18 


E (V/cm) 


<a 

s*W„ 


Figure  6.  Maintenance  Electric  Fields  for  Steady-State 
Microwave  Discharges  Contrasted  with  Breakdown  Values 
(Courtesy  of  S.  C.  Brown1^) 


required  to  initiate  breakdown  without  the  glow  discharge.  The  phenomenological 
interpretation  presented  for  these  results  involves  a  complicated  'precursor' 
mechanism  of  resonance  radiation  that  produces  excited  metastable  and  radiation- 
trapped  states.  It  is  postulated  that  these  excited  states  may  be  ionized  by  elec¬ 
tron  impact  with  electrons  oscillating  in  the  microwave  field  and  that  at  the  level 
where  the  Debye  length  becomes  comparable  to  the  container  dimensions,  a 
charge- separation  mechanism  causes  an  electron  avalanche  to  occur  at  this  re¬ 
duced  breakdown  threshold. 

Bandel  and  MacDonald  have  reported  experimental  measurements  of  the  de¬ 
crease  in  breakdown  field  strength  in  neon  as  a  function  of  increased  preionization 
that  is  the  residual  ionization  from  a  previous  breakdown  pulse.  A  sufficient  in¬ 
terpretation  of  their  results  is  that  the  initial  level  of  ionization  determines  a  cor¬ 
responding  loss  rate  for  electrons  and  that  the  breakdown  field  is  the  field  whose 
ionization  frequency  will  balance  and  infinitesimally  exceed  this  loss  rate.  Since 
the  population  of  metastable  states  is  not  known,  no  conclusions  can  be  drawn  as 
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to  the  mechanism  of  electron-impact  ionization  of  metastables;  however,  the 

mechanism  _*  nnot  be  discounted  as  an  alternative  interpretation  of  the  experiment 

without  additional  experimental  data, 

15 

Recent  studies  of  a  helium  afterglow  irdicate  the  plausibility  of  the  ionizing 
metastable-metastable  collision  as  a  mechanism  to  account  for  discrepancies  in 
experimental  determinations  of  the  electron-ion  recombination  coefficient.  In  a 
model  for  the  afterglow  that  includes  this  mechanism  as  a  source  of  free  electrons, 
the  electron  density  is  directly  proportional  to  the  metastable  density  over  a 
limited  time  interval  in  the  afterglow.  The  congelation  intuitively  expected  be¬ 
tween  density  and  microwave  breakdown  field  is  not  by  itself  strong  evidence  in 

favor  of  any  particular  mechanism. 

3 

Experiments  have  shown  that  microwave  breakdown  occurs  at  a  lower  elec¬ 
tric  field  in  neon-argon  mixtures  than  it  does  in  either  of  these  gases  in  pure 
form.  An  interpretation  i3  that  the  excitation  of  metastable  states  in  the  pure  gas 
constitutes  an  energy  sink  that  requires  additional  energy  input  to  initiate  break¬ 
down.  In  a  Penning  mixture,  this  energy  sink  is  destroyed  by  ionization  of  the 

14 

'impurity'  atom.  Similarly,  in  the  experiment  reported  by  Bethke  and  Ruess, 
it  may  be  hypothesized  that  the  radiative  precursor  destroys  this  energy  sink. 


5.  RESULTS 

The  breakdown  fields  calculated  here  show  a  reduction  in  the  threshold  volt¬ 
age  for  breakdown  as  diffusion  is  modified  by  charge-separation  effects.  In  these 
calculations  the  expressions  involving  the  mobilities  and  diffusion  coefficients  for 
electrons  and  ions  are  replaced  by  constants  A  and  B  so  that  the  transition  is 
forced  to  occur  at  the  same  electron  density.  The  pressure-dependence  of  the 
shape  of  the  transition  is  taken  care  of  by  including  the  factor  1/p  in  the  constant 
B .  The  results  obtained  in  this  program  are  compared  with  the  experimental  re- 
suits  of  Bandel  and  MacDonald  in  Figures  7(a)  and  (b).  At  both  of  the  pressures 
used,  it  will  be  noted  that  the  reduction  in  the  breakdown  field  reaches  a  limit  after 
about  10  electrons/cm  ,  thus  showing  good  agreement  between  experiment  and 
theory.  This  limit  reflects  the  ratio  of  A/B  in  the  calculation  which  is  pressure- 
dependent  through  ,i_  in  (20).  In  the  6.5-Tori  case  the  calculated  curve  closely 

approximates  the  experimental  curve.  For  the  0. 35-Torr  case,  however,  the 

7  8  3 

experimental  results  show  a  sharper  'knee'  between  10  and  10  electrons/cm 
than  the  calculated  values  do.  A  somewhat  higher  electron  temperature  could  be 


15.  Johnson,  A.  W, ,  and  Gerardo,  J.  B.  (1972)  Recombination  and  ionization  in  a 
molecular- ion-dominated  helium  afterglow,  Phys,  Rev.  A5:1410. 
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expected  at  the  lower  pressure  and  hence,  from  (14),  A  and  B  would  have  smaller 
values  as  the  electron  temp,  roture  increased  (the  ratio  of  A  to  B  remains  un¬ 
changed  by  electron  temperature).  As  n  consequence,  a  higher  value  of  electron 
density  in  (16)  would  bo  required  to  cause  :,B  to  become  significant  with  respect 
to  unity  and  hence  the  breakpoint  would  occur  at  higher  electron  Densities,  and 
thus  possibly  more  in  agrt  nment  with  experiment. 
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6.  CONCLUSIONS 


As  a  result  of  these  calculations,  we  can  consider  the  ambipolar  diffusion 
mode  as  a  consistent  explanation  of  the  mechanism  producing  experimental  re¬ 
sults  such  as  those  of  Bandel  and  MacDonald;  however,  further  experimental 
results  are  necessary  to  identify  this  mechanism  as  dominant  in  the  preionization 
case.  The  ratio  of  breakdown  fields  for  the  un-ionized  and  the  preionized  cases 
is  an  observable  that  does  not  require  accurate  absolute  measurement  of  the  break¬ 
down  electric  field  or  the  preionization  electron  density,  and  is  consequently  em¬ 
phasized  here.  An  independent  measurement  of  the  actual  metastable  density  in 
the  preionization  case  is  considered  rifecessary  to  firmly  establish  the  mechanism 
responsible  for  the  phenomenon  of  Bethke  and  Ruess  and  the  general  problem  of 
breakdown  in  the  presence  of  preionization.  Calculations  are  based  on  an  assumed 
electron  temperature  that  was  selected  as  appropriate  at  moderate  pressures. 

The  assumption  of  a  single  value  of  electron  temperature  over  three  decades  of 
pressure  is,  however,  unjustifiable.  Electron  temperature  variations-and  if 

possible,  as  a  function  of  afterglow  time  and  electron  density-need  to  be  examined, 

-2 

especially  since  the  constants  A  and  B  are  roughly  proportional  to  T_  ,  and  the 
transition  onset  is  thus  very  sensitive  to  electron  temperature. 
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